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Abstract 

The production of D*±(2010) mesons in e^p scattering in the range of exchanged 
photon virtuahty 0.05 < < 0.7 GeV^ has been measured with the ZEUS de- 
tector at HERA using an integrated luminosity of 82 pb~^. The decay channels 
D*~^ D^7T~^ with D° — > K~n~^ and corresponding antiparticle decay were used 
to identify D* mesons and the ZEUS beampipe calorimeter was used to identify 
the scattered electron. Differential D* cross sections as functions of Q^, inelastic- 
ity, y, transverse momentum of the D* meson, pt{D*), and pseudorapidity of the 
D* meson, r]{D*), have been measured in the kinematic region 0.02 < y < 0.85, 
1.5 < pt{D*) < 9.0 GeV and \ri{D*)\ < 1.5. The measured differential cross 
sections are in agreement with two different NLO QCD calculations. The cross 
sections are also compared to previous ZEUS measurements in the photoproduc- 
tion and DIS regimes. 
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1 Introduction 



The production of charm quarks at HERA has been studied both in deep inelastic scat- 
tering (DIS) [1-5] and photoproduction [6-10]. In general, reasonable agreement is seen 
with next-to-leading-order (NLO) QCD predictions. 

This paper presents measurements of the D* cross section in the range 0.05 < < 
0.7 GeV^. The beampipe calorimeter of ZEUS [11,12] was used for the measurement of 
the scattered lepton, which allows the first measurements of the transition region between 
photoproduction (photon virtuality, ~ GeV^) and DIS (Q^ > IGeV^). The cross 
sections are compared to the predictions of two different NLO QCD calculations, one de- 
signed for DIS, the other for the photoproduction region. This paper investigates whether 
the calculations remain valid in this transition region. 

2 Experimental set-up 

This analysis was performed with data taken from 1998 to 2000, when HERA collided 
electrons or positrons^ with energy Eg = 27.5 GeV with protons of energy Ep = 920 GeV. 
The combined data sample has an integrated luminosity of £ = 81.9 ± 1.8 pb~^. 

A detailed description of the ZEUS detector can be found elsewhere [13]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [14], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD consists 
of 72 cylindrical drift chamber layers, organized in nine superlayers covering the polar- 
angle^ region 15° < 6* < 164°. The transverse- momentum resolution for full-length tracks 
is a{pT)/pT = 0.0058pr © 0.0065 © 0.0014/pt, with pr in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [15] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a{E)/E = 0.18/^/E for electrons and 
(t{E)/E = 0.35 /Ve for hadrons, with E in GeV. 

^ Hereafter, both electrons and positrons are referred to as electrons. 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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The scattered electron was detected in the beampipe calorimeter (BPC). The BPC allowed 
the detection of low-Q^ events, where the electron is scattered through a small angle. The 
BPC was used in previous measurements of the proton structure function, F2, at low 
[11,12]. It originally consisted of two tungsten-scintillator sampling calorimeters with 
the front faces located at Z = —293.7 cm, the centre at y = 0.0 cm, and the inner edge 
of the active area at X = ±4.4 cm, as close as possible to the electron-beam trajectory. 
At the end of 1997 one of the two BPC calorimeters was removed; hence, for the analysis 
in this paper, only the calorimeter located on the +X side of the beampipe was utilised. 
It had an active area of 12.0 x 12.8 cm^ in X x y and a depth of 24 radiation lengths. 
The relative energy resolution as determined in test-beam measurements with 1-6 GeV 
electrons was ^EjE = 17%/ ^E (GeV). 

The luminosity was measured from the rate of the bremsstrahlung process ep — * e'~fp, 
where the photon was measured in a lead-scintillator calorimeter [16] placed in the HERA 
tunnel at Z=-107m. 

A three-level trigger system was used to select events online [13,17]. At all three levels, 
the event was required to contain a scattered electron candidate in the BPC. Additionally, 
at the third level, a reconstructed D* candidate was required for the event to be kept for 
further analysis. The efficiency of the online D* reconstruction, determined relative to an 
inclusive DIS trigger, was above 95% [5]. 

3 Kinematic reconstruction and event selection 

Deep inelastic electron-proton scattering, ep eX, can be described in terms of two 
kinematic variables, chosen here to be y and Q^, where y is the inelasticity. They are 
defined as = —q^ = —{k — k'Y and y = Q^/{2P ■ q), where k and P are the four- 
momenta of the incoming electron and proton, respectively, and k' is the four-momentum 
of the scattered electron. The inelasticity, which is the fractional energy transferred to the 
proton in its rest frame, is related to the Bjorken scaling variable x and by = sxy, 
where s = AE^Ep is the square of the electron-proton centre-of-mass energy of 318 GeV. 

The values of y and were calculated using the measured electron scattering angle 
and the energy deposited in the BPC as detailed in a previous analysis [11], which also 
describes the method used for the energy calibration of the BPC. A time dependent 
re-calibration of the energy response was necessary [18], as radiation damage of the scin- 
tillator resulted in a degradation of about 10% by the end of the 2000 running period. 

A series of cuts was applied to reject background. The events were required to have a pri- 
mary vertex within 50 cm in Z of the nominal interaction point. The electron candidates 
in the BPC were required to have -Ebpc > 4 GeV, as the trigger efficiency is low below this 



2 



energy. The electron impact point on the face of the BPC was required to be more than 
0.7 cm from the inner edge to ensure good shower containment. Photoproduction events 
were efficiently rejected by requiring the events to have 35 < E — Pz < Q5 GeV, where 
E — Pz = Yliii^ ^ Pz)i is summed over all CAL deposits, including the scattered electron 
candidate in the BPC. Finally, events with an additional well-reconstructed electron can- 
didate in the CAL with energy greater than 5 GeV were rejected to reduce background 
from DIS events with Q"^ > I GeV^. 

The measured kinematic region in y and was restricted to the range of high acceptance, 
0.02 < y < 0.85, 0.05 < Q'^ < 0.7 GeV^. With these cuts, the reconstructed invariant 
mass of the hadronic system, W, lies between 50 and 300 GeV, with a mean of 190 GeV. 

4 Selection of D* candidates 

The D* mesons were identified using the decay channel D*'^ D^T^t with the subsequent 
decay K~tt^ and the corresponding antiparticle decay chain, where 7r+ refers to a 

low-momentum ("slow") pion accompanying the D^. 

Charged tracks measured by the CTD and assigned to the primary event vertex^ were 
selected. The transverse momentum was required to be greater than 0.12 GeV. The px 
cut was raised to 0.25 GeV for a data subsample corresponding to (16.9 ± 0.4) pb~^, for 
which the low-momentum track-reconstruction efficiency was lower due to the operating 
conditions of the CTD [19]. Each track was required to reach at least the third superlayer 
of the CTD. These restrictions ensured that the track acceptance was high and the mo- 
mentum resolution was good. Tracks in the CTD with opposite charges and transverse 
momenta px > 0.45 GeV were combined in pairs to form candidates. The tracks were 
alternately assigned the kaon and the pion mass and the invariant mass of the pair, M^n, 
was determined. Each additional track, with charge opposite to that of the kaon track, 
was assigned the pion mass and combined with the D°-meson candidate to form a D* 
candidate. 

A mass window for the signal region of the varying from 1.82 < Mktt < 1.91 GeV 
to 1.79 < M/^TT < 1.94 GeV was used, reflecting the dependence of the CTD resolution 
on pt{D*). The signal region for the reconstructed mass difference AM = {Mk-k-k^ — 
MKn) was 0.1435 < AM < 0.1475 GeV. The requirement of pt{D*) / E^>^"° > 0.1 was 
also applied, where E^^^^° is the transverse energy outside a cone of ^ = 10° deflned 
with respect to the proton direction. This cut rejects background without signiflcantly 
affecting the signal. 

^ The resolution of such tracks is not good enough to separate primary and secondary vertices from c 
and b hadron decays. 
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The D* mesons were selected in the kinematic region 1.5 < pt{D*) < 9 GeV and \'r]{D*)\ < 
1.5. The AM distribution for events with an electron reconstructed in the BPC is shown 
in Fig. 1. To extract the number of D* mesons, the AM distribution was fit using 
an unbinned likelihood method, with a Gaussian to describe the signal and a threshold 
function to describe the combinatorial background. A first estimate of the background 
was given by D* candidates with wrong-sign combinations, in which both tracks forming 
the candidates have the same charge and the third track has the opposite charge. 
These are shown as the shaded region in Fig. 1. The number of D* mesons obtained from 
the fit was N{D*) = 253 ± 25. 



5 Acceptance corrections and systematic uncertain- 
ties 

The acceptances were calculated using the Herwig 6.1 [20] and Rapgap 2.08 [21] Monte 
Carlo (MC) models. Both models simulate charm and beauty production and include 
contributions from both direct and resolved photoproduction. In direct photoproduction 
the photon participates as a point-like particle in the hard scattering process, while in 
resolved photoproduction a parton in the photon scatters on a parton in the proton. The 
generated events were passed through a full simulation of the detector, using Geant 
3.13 [22] and then processed and selected with the same programs as used for the data. 
The CTEQ5L [23] parton density function (PDF) was used for the proton and GRV- 
LO [24] was used for the photon. The charm-quark mass was set to 1.5 GeV. 

The Herwig predictions are in good agreement with the data distributions for both the 
scattered lepton and hadronic variables and so this Monte Carlo was used to correct the 
data for detector effects. For the kinematic region of the measurement 0.05 < < 
0.7GeV^ 0.02 <y < 0.85, 1.5 < Pt{D*) < 9 GeV, and \v{D*)\ < 1.5 the acceptance was 
(1.11 ± 0.03)%. This includes the geometrical acceptance of the BPC, which was about 
9%, and the reconstruction efficiency for the D* decay chain. 

The Rapgap MC gives a similarly good representation of the data and was used to 
estimate part of the systematic uncertainties, as described below. 

The differential cross section for a given observable Y was determined using 

da _ N 

dY ~ A-C-B-AY' 

where is the number of D* events in a bin of size AY, A is the acceptance (which takes 
into account migrations and efficiencies for that bin) and C is the integrated luminosity. 
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The product, B, of the appropriate branching ratios for the D* and decays was set to 
(2.57 ±0.05)% [25]. 

The systematic uncertainties of the measured cross sections were determined by changing 
in turn the selection cuts or the analysis procedure within their uncertainties and repeating 
the extraction of the cross sections [26]. The major experimental sources of systematic 
uncertainty were (the variation of the total cross section is given in parentheses): the 
BPC alignment (^3^%) and energy scale (^i'2%); the uncertainty in the CTD momentum 
scale (1?:5%) and the CAL energy scale (±1%); the pt{D*) / E^>^^° cut (^?:?%) and the D* 
signal extraction The uncertainty due to the MC model (148%) was determined 

by using Rapgap to evaluate the acceptance correction rather than Herwig, as well as 
by varying the fraction of resolved and direct photoproduction processes in the simulation. 

All the above errors were added in quadrature separately for the positive and negative 
variations to determine the overall systematic uncertainty. The overall normalisation has 
additional uncertainties of 2.2% due to the luminosity measurement and 2.0% due to 
knowledge of branching ratios. These are included in the error quoted for the total cross 
section but not in the systematic uncertainties of the differential cross sections. 

6 Theoretical predictions 

Two different calculations were used to evaluate the theoretical expectation for charm 
production. 

The HVQDIS program [27] implements an NLO calculation of charm production in DIS. 
At low Q^, the hadron-like structure of the photon, not included in HVQDIS, is needed 
to regularise the NLO calculation. Therefore predictions from this program are expected 
to lose accuracy in the limit 0. The ZEUS measurements of D* production in DIS 

for > 1.5 GeV^ are in good agreement with the HVQDIS prediction [5]. 

The FMNR program [28] implements an NLO calculation of charm photoproduction which 
includes the hadron-like component of the photon. Electroproduction cross sections can 
be obtained with FMNR using the Weizsacker-Williams approximation [29] and are there- 
fore expected to be reliable only at low Q^, where this approximation is valid. The FMNR 
predictions are in reasonable agreement with ZEUS measurements of D* photoproduc- 
tion [7], considering the theoretical uncertainties. 

It is therefore interesting to see whether these calculations are able to reproduce the data 
in the transition region between photoproduction and DIS. The following parameters 
were used in the calculations for both programs. They were chosen to be the same as 
in a previous publication [5]. A variant of the ZEUS-S NLO QCD global fit [30] to 
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structure-function data was used as the parameterisation of the proton PDFs. This fit 
was repeated in the fixed-flavour-number scheme, FFNS, in which the PDF has three 
active quark flavours in the proton, and Aq^j-, is set to 0.363 GeV. The mass of the 
charm quark was set to 1.35 GeV. The renormalisation and factorisation scales were set 
to fiji = fip = a/Q^ + 4m^ in HVQDIS, while for FMNR they were set to the usual 
choice of yU/j = /ii? = \/px + 5 where is the average transverse momentum squared 
of the charm quarks. The charm fragmentation to a D* is carried out using the Peterson 
function [31]. The hadronisation fraction, /(c D*), was taken to be 0.238 [32] and the 
Peterson parameter, e, was set to 0.035 [33]. The parameters used here for the FMNR 
calculation are different from those used in a previous photoproduction analysis [7] (which 
used rUc = 1.5 GeV) leading to a 20% larger predicted photoproduction cross section. 

For the FMNR calculation the electroproduction cross section, aep, was obtained from the 
photoproduction cross section, a^p(W), using 



where 



max^ 



271 



i + (i-yr), QLx .^f 1 



ln^-2m,y 



y Qmin vQmin Qmi: 

is the photon flux and ymin, 2/max, QmiiT Qma.x define the measurement range in y and Q^. 

The NLO QCD predictions for D* production are affected by systematic uncertainties, 
which were also evaluated as in a previous ZEUS paper [5]^. The sources of system- 
atic uncertainties on the total cross section are: charm quark mass (^13% for HVQDIS, 
+11% for FMNR); renormalisation and factorisation scale (1}3% for HVQDIS, for 
FMNR); ZEUS PDF (±5%); fragmentation (l6°%). For both programs, the systematic 
uncertainties were added in quadrature and are displayed as a band in the figures. 

Theoretical calculations of the total charm cross section in this range can not be 
compared to the present data since D* are only measured in a limited pt and 77 range. 

7 Cross section measurements 

The total cross section for 0.05 < < 0.7GeV^ 0.02 <y < 0.85, 1.5 < Pt{D*) < 9 GeV 
and |77(-D*)| < 1.5 is: 

a{ep eD*X) = 10.1 ± 1.0(stat.)+J](syst.) ± 0.20(BR) nb. 



For the HVQDIS case, following [5], the minimum value for the scales was set to 2mc. 
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where the first uncertainty is statistical, the second from systematic effects (including the 
luminosity uncertainty) and the third from the uncertainties in the branching ratios. 

The prediction from the HVQDIS program is 8.6!l^} g nb, in agreement with the data, while 
the prediction from FMNR is 8.91^4 nb^, also in good agreement. 

The measured differential D* cross sections as a function of Q^, Pt{D*) and ti{D*) for 
the data are shown in Fig. 2 and given in Table 1. The predictions of the NLO calculations, 
including their uncertainties, are shown as bands. The measured differential cross sections 
are well described over the full measured kinematic region by both calculations. 

This analysis was also compared to previous ZEUS measurements of D* production in 
DIS [5] made in the kinematic region 1.5 < < 1000 GeV^ 0.02 < y < 0.7, 1.5 < 
Pt{D*) < 15GeV and \ri{D*)\ < 1.5. In order to directly compare with the results 
presented there, the cross sections were recalculated in the modified kinematic region 
0.02 < y < 0.7. No correction was made for the different upper cut on pt{D*), as the size 
of the effect is ~ 1%. 

For this modified kinematic region, the differential cross section as a function of is 
presented in Fig. 3 and given in Table 2. The systematic errors were assumed to be the 
same as those in the full y range. Figure 3 also shows the previous ZEUS measurement 
and the HVQDIS prediction. The combination of both measurements shows that the 
slope of da/dQ^ changes with Q^; at high the slope is steeper than at low Q^. The 
NLO calculation describes the measured data well over the full range. 

The D* electroproduction cross sections were converted to ■jp cross sections, a^p, in the 
range 1.5 < pt{D*) < 9 GeV and |?7(D*)| < 1.5 (measured in the laboratory frame) 
using the photon flux from Eq. 1. The cross sections are given for W = 160 GeV, which 
corresponds to y = 0.25, close to the mean y of the measured cross sections. The W 
dependence of a^p was evaluated from the data. The uncertainty of this procedure was 
estimated to be 10%. A comparison of the charm photoproduction cross section [7], this 
measurement and the DIS cross sections [5] is shown in Fig. 4. The numbers are tabulated 
in Table 3. The photoproduction point was corrected for the different kinematic range 
and centre-of-mass energy used here using the FMNR program. As can be seen, the 
present measurements are consistent with the photoproduction cross section. A fit using 
a function of the form a{Q^) = SM"^ /{Q"^ + M^), where 5* is the photoproduction cross 
section at = and is the scale at which the 7p cross section changes from the 
photoproduction value to the DIS 1/Q^ behaviour, gives a good description of the data 
over the whole range with 5 = 823 ± 63 nb and = 13 ± 2 GeVl The value of 
found here for charm production is close to 4m^ [34] and significantly larger than that 
found for inclusive data = 0.52 ± 0.05 GeV^ [12]. 

^ The contribution from the hadron-hke component of the photon is 9%. 
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8 Conclusions 



Charm production has been measured as a function of , y, pt{D*) and ri{D*) in the 
kinematic region 0.05 < < 0.7GeV^ 0.02 <y < 0.85, 1.5 < priD*) < 9.0 GeV and 
|77(D*)| < 1.5. These measurements extend the previous ZEUS measurements in DIS to 
lower Q^. The measured differential cross sections are well described by two different 
NLO QCD calculations: one (FMNR) is designed for the photoproduction region; while 
the other (HVQDIS) is designed for DIS. Both calculations predict similar cross sections in 
the intermediate region measured here, which agree well with the measurements. The 
measurements, converted to 7p cross sections, also agree well with the D* photoproduction 
data. 
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bin 

( GeV^) 


da/dQ'^ Astat 
(nb/GeV^) 


A t 


0.05:0.20 


29.1 


±7.2 


+4.3 
—4.1 


0.20:0.35 


15.0 


±2.4 


+1.5 
-1.4 


0.35:0.50 


10.7 


±2.2 


+1.3 
— 1.1 


0.50:0.70 


7.1 


±2.3 


+1.6 
-0.8 


y bin 


da/dy 


Astat 

(nb) 


A . 

"syst 


0.02:0.15 


34.2 


±6.7 


+7.5 
—7.5 


0.15:0.30 


19.5 


±3.8 


+2.7 
-2.1 


0.30:0.50 


10.7 


±2.1 


+1.2 
— 1.1 


0.50:0.85 


3.8 


±1.1 


+0.8 
—0.8 


v^(D*) bin 
( GeV) 


da/dpT{D*) Astat 

( nb/ GeV) 


A t 

^syst 


1.5:2.5 
2.5:3.8 


6.8 
2.2 


±1.7 
±0.3 


+1.0 
-0.9 
+0.2 
-0.2 


3.8:5.0 


0.53 


±0.11 


+0.02 
-0.02 


5.0:9.0 


0.11 


±0.02 


+0.01 
-0.01 


?7(D*) bin 


da/dr]{D*) 


Astat 

(nb) 


Asyst 


-1.5: -0.5 


3.4 


±0.6 


+0.7 
-0.7 


-0.5: 0.0 


4.1 


±0.9 


+0.5 
-0.4 


0.0: 0.5 


2.9 


±0.8 


+0.3 
-0.3 


0.5: 1.5 


3.3 


±0.7 


+0.4 
-0.3 



Table 1: Measured differential cross sections as a function of Q^, y, pt{D*) 
and r]{D*) for 0.05 < < 0.7 GeV^, 0.02 < y < 0.85, 1.5 < Pt{D*) < 9 GeV 
and |?7(D*)| < 1.5. The statistical and systematic uncertainties are shown sepa- 
rately. The normalisation uncertainties from the luminosity measurement and the 
branching ratios are not included in the systematic uncertainties. 
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g2 bin 


da/dQ"^ Astat 


( GeV') 


( nb/ GeV^) 


0.05:0.20 


30.0 ±7.2 


0.20:0.35 


14.0 ±2.3 


0.35:0.50 


10.3 ±2.1 


0.50:0.70 


6.9 ±2.3 



Table 2: Measured differential cross sections as a function of Q'^ for 0.02 < y < 
0.7, 1.5 < pt{D*) < 9 GeV and \r]{D*)\ < 1.5. The systematic uncertainties are 
assumed to be the same as those for the kinematic range 0.02 < y < 0.85. 
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( GeV') 




^stat 

(nb) 


^syst 


~ 


729 


±46 


+110 
-92 


0.10 


710 


±170 


+200 
-200 


0.26 
0.42 
0.59 


810 
940 
890 


±130 
±200 
±290 


+180 
-180 
+260 
-260 
+370 
-360 


2.7 

7.1 


741 
506 


±31 
±27 


+95 
-100 
+81 
-59 


14 


408 


±22 


+64 
-47 


28 


278 


±13 


+36 
-33 


57 


152 


±13 


+24 
-24 


130 


64 


±9 


+ 14 
-11 


450 


21 


±5 


+6 
-11 



Table 3: 7p cross sections for D* production in the range 1.5 < Pt{D*) < 9 GeV 
and |?7(-D*)| < 1.5 as a function of Q'^ for W = 160 GeV. The values at Q"^ ^ 
and for > 2.7 GeV"^ are obtained from previous photoproduction [7] and DIS 
measurements [5] in the range 1.5 < Pt{D*) < 15 GeV and \'r]{D*)\ < 1.5. 
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Figure 1: The distribution of the mass difference, AM = M^KTrTis) — M{K'k), 
for D*^ candidates with a measured scattered electron in the BPC. The histogram 
shows the AM distribution for wrong charge combinations, normalised to the data 
in the region 0.151 < AM < 0.167. The normalisation factor is 1.07. The solid 
curve is the result of the fit described in the text. 
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Figure 2: Differential D* production cross sections as a function of (a) , (b) y, 
(c) Pt{,D*) and (d) r]{D*) compared to the HVQDIS and FMNR NLO predictions. 
Data are represented by points. The inner error bars are the statistical errors of 
the measurement while the open error bars are the sum of statistical and system- 
atic uncertainties added in quadrature. The shaded area indicates the theoretical 
uncertainties obtained by variation of the HVQDIS parameters. The dashed and 
dotted lines represent the central value of the FMNR calculation and its uncertainty, 
respectively. 
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Figure 3: The D* production cross section as a function of in the kine- 
matic region 0.02 < y < 0.7, 1.5 < Pt{D*) < 9 GeV and |?7(-D*)| < 1.5 for 
this measurement (BPC) and previous results on D* production in DIS [5] (for 
1.5 < pt{D*) < 15 GeV), compared to the HVQDIS NLO prediction. The data 
are represented by points. The inner error bars are statistical while the open error 
bars are the sum of statistical and systematic uncertainties added in quadrature. 
The shaded area indicates the theoretical uncertainties obtained by variations of the 
HVQDIS parameters. 
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Figure 4: The 7p cross section for D*^ production in the range 1.5 < Pt{D*) < 
9 GeV and \ri{D*) \ < 1.5 as a function of from this paper (BPC), compared with 
previous results on D* production in DIS [5] and photoproduction [7] for 1.5 < 
Pt{D*) < 15 GeV and \ri{D*)\ < 1.5. The data are represented by points. The 
inner error bars are statistical while the open error bars are the sum of statistical 
and systematic uncertainties added in quadrature. The photoproduction point is 
drawn at = 0.003 GeV^ for convenience. The curve shows a fit to the data 
described in the text. 
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